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on the slope of the primordial power spectrum n. It is possible to analytically calculate the contribution 
over the full range of scales and redshifts, correctly taking into account fluctuation growth and damping as 
well as thermalization processes. We find that the 95% upper limit is weakly dependent on cosmological 
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for fio < 1 in a flat fio + = 1 universe. 
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1. Introduction 

Two years of data from the COBE DMR experiment has yielded a firm large scale normalization of 
the power spectrum and indications that the spectral index may be n > 1. The temperature fluctuations 
on a 10° scale have been measured to be (AT/T) rms (10°) = 1.12 ± 0.10 x 10" 5 , and the best fit spectral 
index is n = 1.42 ±0.35 (Bennett et al. 1994, Wright et al. 1994b). The Tenerife experiment (Hancock et al. 
1994) measures a high fluctuation level on scales of 5° that, together with the first year COBE DMR result, 
already requires n > 0.9. Furthermore, the COBE FIRAS experiment has stringently constrained spectral 
distortions of the chemical potential type to be |yu| < 3.3 x 10 -4 (95% CL, Mather et al. 1994). We shall 
show below that, regardless of any other constraints from large-scale structure observations on n, a value of 
n > 1.5 results in excessive //-distortions of the spectrum for adiabatic perturbations. 

On small scales, photon diffusion dissipates fluctuations and thereby creates spectral distortions (Sun- 
yaev k Zel'dovich 1970; Daly 1991; Barrow k Coles 1991). Daly (1991) noted that diffusion effectively 
superimposes blackbodies from regions with different temperatures, resulting in a distortion initially of the 
Compton-t/ type with amplitude approximately the square of the temperature dispersion (Zel'dovich, Illari- 
onov, k Sunyaev 1972). Wright et al. (1994a) applied this method to the recent COBE results, but did not 
give a proper treatment of the matter and radiation transfer functions, which relate the fluctuations at any 
given time to the primordial power spectrum. 

Compton scattering processes the initial distortion into a jj, distortion in a manner which depends 
sensitively on the balance between the perturbed energy and number density of photons. Because both 
quantities are of order the variance of the temperature distribution, the calculation is effectively second 
order (Hu, Scott, k Silk 1994). The superposition function must take into account complications such as the 
energy associated with the coherent motions and pressure gradients of the fluctuations. Since the behavior 
of jj, distortions is independent of the detailed form for the initial spectral perturbation, we instead base our 
treatment on an accurate account of the number and energy density alone (Sunyaev k Zel'dovich 1970). 

In this Letter, we improve on the results of Daly (1991) and Barrow k Coles (1991) by explicitly 
including the effect of thermalization and the damping of fluctuations from all wavelengths at all times; 
prior treatments considered the effect to arise from a single wavelength instantaneous perturbation to the 
spectrum at the thermalization epoch. Furthermore, we normalize the power spectrum to the COBE DMR 
detection at 10°. Our treatment uses the correct transfer functions both at small scales, to determine the 
physics of damping, and at large scales, to fix the normalization. We find that the resulting constraint on the 
spectral index n of the matter power spectrum P(k) = Ak n is significantly more stringent than previously 
inferred. 

2. Spectral Distortions 

Before redshift z y = 2.15 x lO 4 ©^ 2 (^l B h 2 )~ 1/2 (Burigana et al. 1991), Compton scattering rapidly 
establishes equilibrium with the electrons, following an injection of energy or creation of photons. Here 
h = iZo/100 km s" 1 Mpc" 1 , © 2 . 7 = T /2.7K, and the present temperature of the CMB is T = 2.726±0.005K 
(Mather et al. 1994). Equilibrium requires that the photon distribution function / is given by a Bose-Einstein 
spectrum, / = [exp(yU + hv/kTe)] -1 , where T e is the temperature of the electrons. After z y , equilibrium 
cannot be established, and any further energy injection will cause a spectral distortion which is approximately 
of the Compton-t/ type. True thermal equilibrium with / as a Planck spectrum (jj, = 0) cannot in general 
be attained since Compton scattering does not change the number of photons in the spectrum. However, 
higher order processes such as bremsstrahlung (e~ + X — > e~ + X + 7, X is any ion) and double Compton 
scattering (e~ + 7 — > e~ + 7 + 7) can create photons and establish thermal equilibrium. For a universe 
with £lsh 2 ;$ 0.1 as required by nucleosynthesis, double Compton scattering dominates the thermalization 
process (Danese k de Zotti 1982, Hu k Silk 1993). 

We can therefore follow the creation and evolution of spectral distortions by balancing the energy 
injection caused by dissipation with the photon creation caused by double Compton scattering. The com- 
putational technique, pioneered by Sunyaev k Zel'dovich and refined by Danese k de Zotti (1977), involves 
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solving the evolution equations for the photon number and energy density. Assuming that photon creation 
is purely due to double Compton scattering, we obtain (Hu & Silk 1993), 



da a Q 

at %dc P~f 

where toe = 2.09 x 10 33 (1 — Y p /2)~ 1 (£lBh 2 )~ 1 &2^ 2 z~ 9 ^ 2 s is the double Compton thermalization time, 
Y p ?s 0.23 is the primordial mass fraction in helium, and Q / p~ t is the fractional rate of energy injection. 
Equation (1) has the simple solution 

p « 1.4 f (Zy) dt^ exp[-(z/z„) 5 / 2 ], (2) 
Jo P~t 

where z^ = 4.09 x 10 5 (1 - Y I) /2)~ 2 l"'<d\[^(Q, B h 2 )~ 2 l"' is the thermalization redshift for double Compton 
scattering. 

3. Damping of Small— Scale Fluctuations 

The damping of primordial density perturbations before recombination injects energy into the CMB 
spectrum and leads to spectral distortions. Furthermore, given the current stringent limit on spectral 
distortions, the arguments of Daly (1991) show that photon diffusion as opposed to nonlinear dissipation is 
the dominant mechanism for the creation of distortions. 

Let us now examine the photon diffusion mechanism. At early times, photons and baryons are tightly 
coupled and behave as a single viscous fluid. Any fluctuations in the photon-baryon fluid, whether adiabatic 
or isocurvature, will oscillate as a sound wave inside the photon-baryon Jeans length. Solving the Boltzmann 
equation for the photons in the tight coupling limit leads to for the photon energy density perturbation 
(Peebles & Yu 1970), 

A(k, t) = Aj(fc) cos(w^) exp(-/dr ? r), (3) 

where Aj(fc) is the amplitude of the fc-mode upon entering the Jeans length, r\ = J dt/ a is the conformal 
time, a = 1/(1 + z) is the scale factor, the dispersion relation for acoustic waves is to = c s k, and (Peebles 
1980) 

r(„) = p * 3 +w+w\ (4) 

[l> 6n e a T a (R+l) 2 1 ' 

Here R = Zp},/Ap- n n e is the electron number density, a? is the Thomson cross section, and c = 1. This 
damping of the acoustic waves will dump energy into the microwave background. The average energy density 
in a plane acoustic wave is given by p s ?a Ay&cJ (S 2 b ) , where pjj, = p 7 + pj, and Sjj, are the density and density 
perturbation in the photon-baryon fluid, and the brackets denote an average over an oscillation of the acoustic 
wave. Since p distortions arise at z > z y > z eq (= 2.50 x 10 4 fi /j 2 ©2.7)> 

we can take the radiation-dominated 

limit of these equations, R = 0, c 2 s = 1/3, and 

(S 2 b ) « <A 2 (M)> = \A 2 (k)ex P [-(k/k D (t)) 2 ], (5) 

where k D (z) = 2.34 x 1O- 5 2 .7(1 - Y p /2) 1 ' 2 (^l B h 2 ) 1 ' 2 z 3 ' 2 Mpc" 1 . Therefore, the rate of fractional energy 
injection 

Q(t) = v l d(A 2 (fc,t)) 
Pi £ 3 dt 

Inserting this expression into equation (2) allows us to calculate the spectral distortion from photon diffusion. 
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4. Constraining the Power Spectrum 



Hitherto, our arguments have been fully general and apply to any fio or spectrum of fluctuations. We 
will now specialize the argument to the case of primordial adiabatic fluctuations with P(k) oc k n in a flat 
fio + f^A = 1 universe. Adiabatic perturbations in the photon-baryon fluid grow as a 2 outside the Jeans scale, 
which is nearly equal to the horizon scale in the radiation-dominated epoch. The time at which a given mode 
k crosses the horizon is an oc Therefore, in addition to the primordial spectrum, the fluctuations gain 

a factor k~ 2 . This specifies a transfer function such that A 2 (k,t) = T 2 (k,t)Ak n , which, including photon 
diffusion within the Jeans scale in the radiation-dominated epoch, is given by (Kodama & Sasaki 1986), 

T(k, t) = -5a eq k 2 eq k~ 2 cos(w^) exp[-(k/k D ) 2 /2], (7) 

where k eq = (Ha) eq is the comoving wavenumber that crosses the horizon at matter-radiation equality 
a eq = 1/(1 + z eq ). The transfer function T(k,t) has been normalized so that at present on large scales, 
P(k) = Ak n . Therefore for adiabatic fluctuations, Aj is defined as 

A 2 (k) = 2ha 2 eq ki q Ak n -\ (8) 

The normalization constant A can be determined by the COBE detection at 10°. 

We now make the further assumption that the COBE DMR detection probes potential fluctuations on 
the last scattering surface via the Sachs-Wolfe effect (Sachs & Wolfe 1967), assumed to arise from scalar 
perturbations alone. We can then write 

/AT\ 2 1 °° 

("T ) = ^£^+WQ, W 

V 7 1=2 

where the COBE window function is Wi = exp [-£(£+ l)cr 2 ], with a = 0.0742 being the gaussian width of 
the 10° FWHM beam and (Efstathiou, Bond k White 1992) 

Ct = V 1 )a i — t — ^-4av x , (10) 
8^F^- 1 r(^)r(£+^) *' v ; 

where V x is the volume over which the universe is assumed to be periodic. Equation (10) takes the radiation 
transfer function at large scales fully into account for an fio = 1 universe. For fio < 1 with a cosmological 
constant, there is a small correction due to the integrated Sachs- Wolfe effect which boosts the low multipoles. 
Given the present status of the observational uncertainties, its inclusion will not significantly affect our 
constraints (Efstathiou, Bond, k. White 1992). 

We are now in a position to evaluate equation (2) explicitly. Both the sum over k and the integral over 
time can be performed analytically for power-law initial spectra leading to, 

^=lAF(n)^[k 3 D A 2 (k D )] z ^ (11) 

where 

F(n) = ±T[(n + l)/2] T[3(n - l)/5, (z y /z,) 5 / 2 ]. (12) 

If n is significantly greater than unity, the incomplete gamma function T(m, x) — > T(m) since Zy/z^ <C 1 and 
F(n) is roughly of order unity. 

It is easy to interpret this result. If n > 1, the smallest waves carry the most energy, and the distortion 
comes almost entirely from the waves that damped at the thermalization epoch. Prior to thermalization, no 
distortion survives due to the rapidity of the double Compton process. In this way, we have not only justified 
the approximations of Daly (1991) and Barrow & Coles (1991) which involve instantaneous injection from a 
single fc-mode, but have rigorously derived the additional n-dependent factor F(n). 



4 



On the other hand if n < 1, the fractional energy injection from dissipation will be a maximum at the 
latest relevant time, i.e. recombination. This implies that the constraint from spectral distortions will come 
from the upper limit on Compton-t/ distortions. Nevertheless, this constraint will not be very powerful since 
the normalization is fixed at large scales. For n < 1, the small-scale power decreases with respect to that 
on large scales, yielding no constraint on a COBE-normalized spectrum. 

Note that the dependence on the cosmological parameters fio, h is quite weak: approximately 

(< 10 7 J 6/5 ) 1 -"4 2 -" )/2 . Hence for n 1, is completely independent of h and f2#. Moreover, yU 
is nearly independent of for all n, since raising makes both the damping length shorter and the 
thermalization redshift smaller. It is also useful to provide an approximate inversion of equation (11): 

ln[g 1 i] - M V/(Ar/T)j .] 
~ + ln[A' 2 (fi B 7 l 2)-i/io(fi 7 l 2)-i/2/(fi )] ^ > 

where we find the constants K\ = 5.6 x 10 -3 and K2 = 8.9 x 10 5 and the small logarithmic correction 
J(f2o) ~ 1 — 0.085 In fio- One can verify that this is an excellent approximation within the range 1.0 < n < 2.0 
and the allowable cosmological parameters. Note that the dependence of n — 1 on yU and the normalization 
is only logarithmic, and its dependence on the cosmological parameters is almost entirely negligible. Even 
relatively large changes in yU or the normalization will not greatly affect the constraint on n. 




FIG 1. Spectral distortion yU, predicted from the full equation (11), as a function of the power index n 
for a normalization at the mean of the COBE DMR detection (AT/T) 10 ° = 1.12 x 10" 5 . With the 
uncertainties on both the DMR and FIRAS measurements, the 95% upper limit is effectively jj, < 
0.63 x 10 -4 (see text). The corresponding constraint on n is quite weakly dependent on cosmological 
parameters: n < 1.54 (h = 0.5) and n < 1.56 (h = 1.0) for fio = 1 and quite similar for 0.2 < fio = 
1 — < 1 universes. These limits are nearly independent of We have also plotted (dashed lines) 
the 95% upper limit on |yu| < 3.3 x 10" 4 for comparison. Note that there is nothing awry with the 
placement of this figure. 

The best fit value of jj, to the spectral data from the FIRAS experiment is jj, = — 1.2 ± 1.1 x 10" 4 (68% 
CL) (Mather et al. 1994). We can actually make a stronger statement than the quoted FIRAS 95% confidence 
limit of |yu| < 3.3 x 10 -4 , by noting that photon diffusion never creates a negative chemical potential and the 
current constraint on positive jj, is somewhat tighter than that on |yu|. Taking into account the COBE DMR 
measurement errors and adopting a 4//K cosmic variance, (AT/T) rms (10°) = 1.12 ± 0.18 x 10 -5 (Bennett 
et al. 1994), we find 



(AT/T) 



, < 5.0 x 10 5 (95%CL). (14) 

10° 
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This would be equivalent to an upper limit of < 6.3 x 10 -5 for a fixed normalization at the mean value 
of the DMR detection. Using (14), we set a limit on the slope n < 1.54 for h = 0.5 and n < 1.56 for 
h = 1.0 (see Fig. 1) for fio = 1 with similar but slightly less stringent limits for cosmological constant- 
dominated universes (fio < 1)- These constraints are nearly independent of and do not change within 
the nucleosynthesis bounds of 0.011 < 9, B h 2 < 0.016 (Walker et al. 1991, Smith et al. 1993). 

5. Discussion 

If the present indications of a large n at COBE DMR scales are borne out by future measurements, 
we should consider how the constraint might be avoided or at least weakened. If the universe is open 
(flo + < 1)) the radiation transfer function at large scales (see equation (10)) responsible for the COBE 
DMR normalization will differ. Significant temperature fluctuations also arise in this model from the inte- 
grated Sachs- Wolfe effect (Sugiyama & Gouda 1992, Hu & Sugiyama 1993) which enhances the large scale 
fluctuations. However, for open universes, there is the more significant effect of an effective cutoff to the 
temperature perturbations at large scale due to the curvature (Wilson 1983; Kamionkowski & Spergel 1993). 
This can naturally lead to an apparent n > 1 slope at large scales without affecting the slope at small scales. 
However, note that the dependence of the constraint on the normalization is only logarithmic. In these 
cases, the upper limit on n, the primordial power index, is roughly the same, but does not correspond to the 
effective n that the COBE DMR experiment measures. 

To change the limit on the primordial power index n, one must change the radiation transfer function 
at small scales where the physics of damping occurs. If the perturbations are not adiabatic, the relation 
given in equation (7) is not valid and the constraint on n will be altered, e.g. for isocurvature fluctuations 
(Daly 1991). These considerations point to a general way of evading the constraints. We have assumed 
that the primordial power spectrum is a pure power law from the current horizon to the damping scale at 
thermalization. If we alter the initial conditions so that the effective n decreases sufficiently with increasing 
k, we can always escape constraints from spectral distortions for any given n at large scales. 

Aside from these considerations, our constraint is difficult to evade since it relies only on the relatively 
well understood physics of photon-baryon systems. Similar limits can be applied from other constraints, 
e.g., primordial black holes (Carr & Lidsey 1993, Novikov et al. 1979) or from avoiding excessive anti-bias 
in galaxies {e.g. Lyth & Liddle 1994). However, such limits could be avoided by appealing to our ignorance 
of the relevant physical processes. By contrast, we believe that the jj, constraint of roughly n < 1.55 is a 
robust limit on pure power laws with adiabatic initial conditions. 

This research has been supported at Berkeley in part by grants from NASA and NSF. 
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FIGURE CAPTIONS 



FIG 1. Spectral distortion predicted from the full equation (11), as a function of the power index n 
for a normalization at the mean of the COBE DMR detection (AT/T) 10 ° = 1.12 x 10" 5 . With the 
uncertainties on both the DMR and FIRAS measurements, the 95% upper limit is effectively jj, < 
0.63 x 10 -4 (see text). The corresponding constraint on n is quite weakly dependent on cosmological 
parameters: n < 1.54 (h = 0.5) and n < 1.56 (h = 1.0) for fio = 1 and quite similar for 0.2 < fio = 
1 — < 1 universes. These limits are nearly independent of We have also plotted (dashed lines) 
the 95% upper limit on \fi\ < 3.3 x 10 -4 for comparison. 
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